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| Abstract. A bipartite quantum state is tomographically faithful when it can be 

çvq ' used as an input of a quantum operation on one of the two quantum systems, such 

ÇN| ■ that the joint output state carries a complete information about the operation itself. 

Tomographically faithful states are a necessary ingredient for tomography of quantum 
operations and for complete quantum calibration of measuring apparatuses. In this 
ÇS| , paper we provide a complete classification of such states for continuous variables in 

C^l ' terms of the Wigner function of the state. For two-mode Gaussian states faithfulness 

i simply resorts to correlation between the modes. 

o 
o 



PACS numbers: 03.65.Wj 
1. Introduction 



Quantum operations describe any kind of physical process affecting quantum states, 



including unitary evolutions of closed systems and non unitary transformations of open 
quantum systems, such as systems interacting with a reservoir, or subjected to noise or 
measurements of any kind. The problem of determining experiment ally the quantum 
operation occurs in different scenarios, typically for quantum calibration of controlled 
transformations PQ and of measuring apparatuses [2] . 

In a naive process tomography one varies the input state over a suitably complete 
set in order to recover enough information about the quantum operation. The tensor 
structure of the bipartite quantum system, however, allows to use a single fixed bipartite 
state that scans the complete set of single-system states in a quantum paral·lel fashion 
PP. The bipartite states that can be used in this way in order to carry a complete 
information of the process are called tomographically faithful [3] . 

The problem of evaluating the faithfulness of a state can be expressed in terms 
of an invertibility condition of a map associated to the state. In many situations it 
is not simple to check such a condition. In this paper we address the continuous 
variables case (i.e. quantum harmònic oscillators), and solve the problem of the complete 
classification of faithfulness in terms of the Wigner function of the state. The use 
of Wigner functions has proved very useful as a generalized phase-space technique to 
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express density operators in terms of c-number functions, thus leading to a considerable 
simplification of the evaluation of quantum dynamics and of expectation vàlues jUE]. 

In this paper we present a general result that provides a necessary and sufhcient 
condition for the faithfulness in terms of Wigner function. Such a condition, as we will 
show, makes use of highly irregular functions, as the customary P-functions in quantum 
òptics. We then specialize our results to the case of Gaussian states of two modes of the 
electromagnètic field. The class of Gaussian states constitutes a fortunate framework 
both for theoreticians and experimentalists, since, on one side all calculations can be 
done analytically, whereas, on the other side, these states are easily generated in a lab, 
using làsers, linear òptics, and parametric amplifiers. We will show that for Gaussian 
states the condition of faithfulness is just the existence of correlations between the two 
modes. 

The paper is organized as follows. In Sec. II we briefly recover the general result 
about the faithfulness of a quantum state, and recali the problem of inversion of a 
special operator associated to the state. We then restate the problem in terms of Wigner 
function of the state, and write a necessary and sufficient condition. The section presents 
some examples of faithful (both entangled and separable) and unfaithful states. In Sec. 
III we simplify the result of Sec. II for the case of two-mode Gaussian states. The 
condition of faithfulness then simply restates as the existence of correlations between 
the two modes. We conclude the paper in Sec. IV, with a summary of results and 
some remarks about the statistical errors that affect the reconstruction of a quantum 
operation, and the connection with quantum images. 



2. Faithfulness in terms of Wigner function 

In mathematical terms, a quantum operation S is described by a completely positive 
map [B]. This can be written in the Kraus form 

£{p) = Y,K nP Kl (1) 

n 

where K n are operators on the Hilbert space 7í of the quantum system. For simplicity 
we will consider quantum operations with the same input and output space TC, and 
that are trace-preserving — the so-called channels — corresponding to the completeness 
relation J2 n Kn^-n = I- The concept of tomographically faithful state |3j relies on using 
a bipartite state R on Ti ® H, such that the output state 

Rs = (S ® I)R (2) 

is in one-to-one correspondence to the quantum operation £. In [3J it was proved that 
a state R is faithful iff the following operator on Ti <8> TC 

R = (ER) T2 E = (R T2 E) T1 (3) 

is invertible. In Eq. Q E = ^\ . denotes the swap operator, and O n denotes 

the partial transposition of the operator O on the /th Hilbert space, / = 1, 2. 
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Using the notation of Ref. [7j for bipartite vectors 

\A)) = y^(n|A|m) |n) ® |m), (4) 

one can generally write a bipartite state in the form 

R = Y J \ A *))((B 1 \- (5) 

From the identity 

A®B\C)) = \ACB T )) (6) 

it follows that 

R = J2AJ®B], (7) 

where the transpose r is defined on the basis chosen in the decomposition of Eq. (jlj). 
Similarly, for a state R written as 

R = ^Ai®B 3 , (8) 

using Eq. Q, one easily writes i? as follows 

A=Ç|S i »(K|, (9) 

where O* denotes the complex conjugation on the fixed basis. Notice that the evaluation 
the operator R does not need the diagonalization of R. Moreover, all the previous sums 
can be suitably replaced with integrals. 

A state R is faithful iff the associated operator R is invertible. We are interested in 
finding the conditions of faithfulness in terms of the Wigner function of the state. For 
simplicity, we consider bipartite states that correspond to two-mode states. However, 
our results are easily generalized to the case of two-party multimode states. 

We recali the Cahill-Glauber formulas [H] between single-mode density matrix p 
and Wigner function W(a, a*) 

W(a,a*) = -Tr [pD (2a) (-lf% (10) 

7T 

p = 2 f d 2 aW(a,a*)D(2a)(-lf a , (11) 



c 



where a* denotes the complex conjugate of a, d 2 a = dR,e(a) rflm(cr), and D(a) = 
e aa^-a*a re p resen ts the displacement operator for the mode a, with [a,a^] = 1. 

For a two-mode bipartite state R we write the Wigner function as a function of two 
complex variables a and f3 by direct generalization of Eq. (fTTj) as follows 



R = A I d l a I d 2 (3W(a,a*,p,p*)D(2a)(-lf a ® D(2(3)(-lf\ (12) 
./ : ./ :• 

where a is a shorthand notation for the annihilation operator otherwise denoted a ® /, 
as well as b for I ® b. In the following we will write the Wigner function shortly 
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as W(a,P) = W(a,a*, P, P*), omitting the dependence on the complex conjugated 
variables. According to Eqs. (JHJ) and (jÜJ), the condition of faithfulness of the state in 
fil 2(1 corresponds to the condition of invertibility of the operator 

R = A I ' d 2 a I d 2 pW(a,p)\D(2p)(-lf a )){(D(2a*)(-lf a \ (13) 
Jc Jc 

= 4(/®(-l) fetfe ) I d 2 a [ d 2 p\V(a,P)\ D(2P) ))((D(2a*)\(I ® (-lf b ) , 
Jc Jc 

where in the second line we used identity (jü|) . Since the set {| D(a)))} is an orthonormal 
basis (in the Dirac sense) for H®H, namely / c D(a)))((D(a) \ = the condition 
of faithfulness in terms of Wigner function is then the following: 



Necessary and sufncient condition for faithfulness: a bipartite state with 
Wigner function W(a,P) is faithful iff one can find a function f(P,j) such that 

d 2 pW(a,P)f(P n ) = 5^ 2 \a- 1 ), (14) 

where ó^^cr) = J c ^ e ^o-*-\*a denotes the Dirac delta over the complex plane. 

Equation (|T4^) should be read in distributional sense. When such a condition is 
satisfied one has 

tfb\ / j2ï I j2., //r .,\ i n/o.,*\\\//n/ojr\ i/r<x c -\\tfb\ 



RT 1 = -{I ® {-ïf b ) / d 2 S / d 2 1 f(6,j)\D(2 1 *)))((D(2S)\(I ® (-1D . (15) 

Jc JC 

2.1. Example 1: twin-beam state 

Consider the twin-beam state that can be easily generated by nondegenerate optical 
parametric amplifiers 

R = (i - A 2 ) | A ata )) ((A ata | , < A < 1 , (16) 

where the parameter A is simply related to the total number of photons n = 2X 2 / (1 — A 2 ). 
The corresponding Wigner function is given by 

W R (a,P) = 4<yl ~ A ^ Ti[\ a ' a D(2a)\ a ' a D(-2p*)] . (17) 

7T 

By normal ordering Eq. f)17j) and a lengthy calculation, one obtains 



4 

W R (a,p) = -^exp 



7T 2 



2(3 : X \ a \ 2 +\P\ 2 ) + -^-(aP + a*P*) 



1 — A 2 ' 1-A 2 

Using the solution (|4*T|) of identity (|T4^) derived in the appendix, the function 7) 
can formally be written as 

mi) = -i^— e^'^e 21 ^ 1 ^ f d^ e \e e ^-^(^-P) (19) 
tt(1 - A 2 ) 2 Jc 

which should be treated as a distribution, in the sense that the integral in £ has to be 
performed after the integration on P of Eq. (|14j) . 
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Notice that in this simple example, the faithfulness is more easily checked using 
Eq. (J2J) and writing immediately 

2.2. Example 2: classically correlated coherent states 

A mixture of correlated coherent states can be easily generated by splitting thermal 
radiation in a 50/50 beam splitter. For such a kind of states we can write 

R= í ^ e -^| 7 >< 7 r, (21) 

where the variance o is related to the total number of photons by n = a/2. The 
corresponding Wigner function is given by 

|2 , 4a 2 



tt^,l) = ^ 2(1 + 2(72) exp 



M + W) + 1 — si( a ^ + a *£) 



(22) 



1 + 2a 2VI 1 y 1 + 2a 2 

Using again the solution (jü]) of identity (fTIJ) derived in the appendix, the function 
7) can formally be written as 

7r(l + 2eHj J c 
and thus the state ({21} is an example of separable faithful state. 

2.3. Example 3: product states 

Consider a product state 

R = p ® d. (24) 

The Wigner function is given by the product of the independent Wigner functions for 
p and a 

W R (a,(3) = W p (a,a*)W a (P,(3*). (25) 
Of course the state R is not faithful, and in fact the condition (fT4*|) can never be satisfied. 

2-4- Example 4- classical correlation between orthogonal states 
Consider the state 

00 

R= (l-A)^À n |n)(n|® 2 , (26) 

n=0 

where \n) denotes the Fock state. From the relation [S] 

(n\D(a)\n) = e-^L n (\a\ 2 ), (27) 
and the identity [Ü] 

f; \«L n (x)L n (y) = T^T ^ {X+V) J o (2^) , (28) 

n=0 ^ ' 
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where L n (x) and Io(x) denote the n-th order Laguerre polynomials and the 0-order 
modified Bessel function, one obtains the Wigner function 

W R (a,(3) = íe-WW/o (^m) ■ ( 2 9) 

Condition (fT4"j) can never be satisfied, since there is no dependence of the Wigner function 
on the phase of (3. In fact, R = R is clearly not invertible, whence the state R is not 
faithful. 



3. Simplificat ion for Gaussian states 

Unfortunately, it is often difhcult to inspect condition (|T4"j) , since it holds more generally 
in a distribution sense. For Gaussian bipartite states, however, it is possible to derive a 
more practical condition in terms of the correlation matrix. 

According to the derivation in the appendix, the term of the Wigner function of a 
Gaussian bipartite state that is relevant for the condition ()14|) is the factor of the form 

g(a, (3) = exp[(Aap + Ba(3*) + h.c] , (30) 

In fact, as shown in the appendix, the condition (|T4j) can be satisfied iff 

|A| 2 -|B|V0. (31) 

In order to clarify the meaning of condition (jUlj) . it is useful to consider the state R in 
terms of the characteristic function T(a,(3) = Tr[RD(a) (g> D(j3)], that corresponds to 
the Fourier transform of the Wigner function, and hence will be Gaussian as well. One 
has 

R= í — í —T(a,{3)D\a)®D\f3). (32) 
Jc n Jc 71 

The operator R then can be written as 

R= — ^r(a,/3)|L>t(/3)))((L>t( a *)|. (33) 
Jc n Jc ^ 

Similarly to Eq. JHJ), R is invertible iff one can find a function 7) such that 

/ d 2 f3T(a } f3)f((3 }7 ) = 5^(a- 1 ) } (34) 
Jc 

and .R -1 can be written as 

R-i= í d 2 5 í d 2 1 f(5, 1 )\D\ 1 *))){{D\5)\. (35) 
Jc Jc 

The same consideration before Eq. (}30|) applies here. The condition (|34j) can be satisfied 
iff |v4| 2 — \B\ 2 ^ 0, where A and B are the coefficients in front of the variables a/3 and 
a(5* of the characteristic function. Using the identities 

A = d 2 al3 T{a ) (3)\ Cí= p = Q - d a T(a, i3)\ a=/3=0 dpT(a, i3)\ a=/3=0 

A* = 9 2 , ( g*r(a,/5)| a=( g =0 - <9 a .r(a,/3)| a=( g = o9 /3 *r(a,/3)| a=( g = o 

B = d 2 a pS(a,(3)\ a=f i = o - d a T(a, P)\ a= /3=odf3*T(a, P)\ a= i3= 

B* = d 2 a , p T(a,(3)\ a= p = o - d a *T(a, (3)\ a= p =0 dpr(a, P)\ a=/3=0 , (36) 
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the condition of faithfulness can be restated in terms of the correlation matrix as follows: 

Necessary and sufncient condition for faithfulness of Gaussian states: a 

bipartite Gaussian state is faithful iff the following condition on the correlations is 
satisfied 

X (R) = (AaW) (Aab) + (Aa^b) (Aa6 f ) ^ 0, (37) 

where for any two operators P and Q 

(APQ) = (PQ)-(P)(Q). (38) 

In terms of the quadratures X c = (c + and Y c = (c — d)/(2i) of the modes 
c = a,b, the correlation x{R) can De rewritten as 

X (R) = \ {(AX a X b ) 2 + (AY a Y b ) 2 + (AX a 7 è ) 2 + (AY a X b ) 2 ) . (39) 

In the examples 1 — 4 given in the previous section, one has x(R) = (ï~\2)2 > °" 4 > > anc ^ 
0, respectively (notice, however, that the state in example 4 is not Gaussian). Using Eq. 
(|39|) the condition x(-R) 7^ shows that bipartite Gaussian states are always faithful, 
except when they are product states (x{R) is the sum of nonnegative terms which all 
vanish when there is no correlation between the two modes). A rigorous relation bewteen 
statistical errors that affect tomographic reconstructions and the strength of correlations 
is in order, but beyond the scope of this paper. 

4. Conclusions 

Tomographically faithful states are a necessary ingredient for tomography of quantum 
operations and for complete quantum calibration of measuring apparatuses. In this 
paper we have provided a complete classification of two-mode faithful states in terms of 
the Wigner function of the state. This classification has been derived from the general 
faithfulness condition resorting to the invertibility of a special operator associated to 
the state. Some examples of faithful states have been presented, both entangled and 
separable, along with examples of not faithful states. For two-mode Gaussian states we 
have shown that faithfulness is simply equivalent to nonvanishing correlations between 
the modes. 

We conclude by noticing that the actual statistical efliciency of a faithful state 
in the tomography of a quantum operation in infinite dimensions is connected to the 
increase of the singular vàlues of the unbounded operator R~ l . Such unboundedness is 
responsible of the increasingly large statistical errors in the Fock representat ion of the 
quantum operation, accounting for the finite experimental data sample used to infer 
information on a infinite set of matrix elements of the quantum operation. As a rule of 
thumb, the statistical efficiency increases for greater correlation x(-R). 

Finally, it is worth mentioning that the framework of quantum images ^0] bears 
a strict analogy with that of the quantum tomography of a channel using an input 
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bipartite state, with the role of the channel here played by the density contour of the 
image analyzed by one of the twin beams from parametric downconversion of vacuum. 
Clearly, when the state is faithful one has quantum imaging on the other beam, and 
our result is consistent with the recent demonstration ^ü] that entanglement is not 
necessary for quantum imaging. In particular, the thermal state split by a beam splitter 
in Eq. (}2T]) is suitable for quantum imaging. 
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Appendix 

We show that for a function of the form 

g(a, f3) = h{a) k{(3) e Aa ^ A * a ^* e B«P*+B*a*p ^ (40) 

with \A\ 7^ \B\, and both h and k (generally not analitycal) invertible functions, the 
following function 

f(P, 7 ) = (|A|2 ~ |jB|2)2 k^W) h^h) e -(l^l 2 +|s| 2 )|7l 2 -(^7 2 +^B* 7 -) e - m » 

j cf £ e (\ A \ 2 +\B\ 2 M\ 2 +(A*B*e+ABÇ*2) e (| J 4| 2 -|B| 2 )($7-í*7·'·) e /3(AC+B·t)-p·(A·t+B?) ^ ^ 



X 

satisfies the identity 

d 2 Pg(a,(3)f((3, 1 ) = 5^(a- 1 ). (42) 



The function /(/?, 7) should be treated as a distribution, in the sense that the integral 
in £ has to be performed after the integration on (3 of Eq. (j42)l . 
One has 



7) = (|A|2 ~ |jB|2)2 h(a) h^h) e-(\A\^\B\>)h\>-(AB^ + A*B*r 2 ) 
X / d 2 Í e (\A\ 2 +\B\ 2 M\ 2 +(A*B*e+ABC 2 ) e (|A| 2 -|Sp)( Í7 -e 7 *) 



d 2 Pg(a,(3)f((3, „ 
c ^ 



(\A\ 2 -\B\ 2 ) 2 



TY 2 



h{a) / i (- 1 )( 7 ) e^l 2+ l B l 2 )(H 2 -l7 



2VU|2_U,|2) +AB ( a 2_ 7 2) +A * B *( a *2_ 7 *2) 



x j d 2 ^ e (|A|2 " |B|2)K(7 - Q) " r(7 *~ Q * )] 

= 6®(a- 1 ), (43) 
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where the integral in d 2 f3 has been performed by using the identity 
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